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Abstract 

In the case of Xe+Sn at 32 MeV/nucleon multifragmentation reaction break-up fragments are 
built-up from the experimentally detected ones using evaluations of light particle evaporation mul- 
tiplicities which thus settle fragment internal excitation. Preeze-out characteristics are extracted 
from experimental kinetic energy spectra under the assumption of full decoupling between frag- 
ment formation and energy dissipated in different degrees of freedom. Thermal kinetic energy is 
determined uniquely while for freeze-out volume - collective energy a multiple solution is obtained. 
Coherence between the solutions of the break-up restoration algorithm and the predictions of a 
multifragmentation model with identical definition of primary fragments is regarded as a way to 
select the true value. The broad kinetic energy spectrum of ^He is consistent with break-up genesis 
of this isotope. 

PACS numbers: 25.70.Pq Multifragment emission and correlations, 24. 10. Pa Thermal and statistical models 
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I. INTRODUCTION 



The possibility to investigate the phase diagram of highly excited nuclear systems made multi- 
fragmentation one of the most lively debated subjects in nowadays physics. However, the approach 
of the thermodynamically relevant break-up stage of the process is difficult from both concep- 
tual and practical perspectives. From the theoretical point of view, the break-up definition of an 
open system requires reconsideration of a fundamental thermodynamical quantity, the volume 
From the experimental point of view, all measured observables are perturbed by non-equilibrium 
phenomena and subsequent disintegration of primary excited fragments. 

Adopting the freeze-out equilibrium hypothesis, the purpose of the present work is to investigate 
to what extent the break-up stage may be restored from limited but relevant experimental data. 
The high quality of the INDRA detector and the wealth of already existing analyses [3, 0, Q] that 
indicate the formation of an equilibrated source make the Xe+Sn at 32 MeV/nucleon multifrag- 
mentation reaction the most suitable candidate for such a study. The present paper follows the line 
opened by Piantelli et al. in Ref. 5| whose goal was that of estimating freeze-out properties such 
as the volume in a fully consistent and model-independent way. Quite remarkably, if constrained 
by all available experimental information, the requirement of optimal fitting of both measured light 
charged particle and fragment kinetic spectra resulted in rather high average excitation energy of 
primary fragments (3.9 MeV per nucleon). While both stick to the break-up concept, the present 
work distinguishes from Ref. [5] by considering that break-up partitions can be traced-back using 
an evaluation of evaporated particle multiplicities jj] and treating all energetic degrees of freedom 
as fully decoupled. 

After a brief review of experimental information, Section II describes the algorithm designed 
to determine break-up fragment partition and average excitation energy. Attempts to identify the 
freeze-out volume {V) and kinetic properties (thermal (K^fi) and collective {Ecoii) energies and flow 
profile acou) of the break-up stage are presented in Section III. Section IV discusses the uniqueness 
of the solution and coherence with results of a multifragmentation model with identical fragment 
definition is thought as a way to select the "real" value and to check the consistency of the results. 
In Section V we review the steps-forward accomplished by different studies meant to characterize 
the break-up stage of the Xe+Sn at 32 MeV/nucleon multifragmentation reaction, their limitations 
and the relevance of present results. 



2 



II. DETERMINING BREAK-UP FRAGMENT PARTITION AND EXCITATION 
ENERGY 



Experimental data files corresponding to Xe+Sn at 32 MeV/nucleon multifragmentation reac- 
tion [j, |3|] are used to reconstruct the thermalized source formed in central collisions. As a first 
step, quasi-complete events where total detected charge exceeds 80% of the summed up charge of 
the target and projectile are selected and the emission isotropy for fragments is ensured by large 
values of the fragment flow angle, 0//oto > 60°. By correcting for detection efficiency, the cross 
section of these selected isotropic central collisions is estimated at 115 it 20 mb 41. Mass of the 
isotopically un-resolved products with Z > 5 is calculated assuming an EAL dependence q]. Then, 
light charged particles {1 < Z < A) emitted at forward {9cm ^ 60") and backward {0cm ^ 120") 
angles in the center of mass of the reaction are rejected as partly created at pre-equilibrium. There- 
fore, light charged particles emitted outside this rejection zone are doubled under the assumption 
of isotropic emission. Neutrons (which are undetected) are added to each event to conserve the 
N/Z of the total system. At this stage the source characterized by < ^4 >=198 and < Z >=83 is 
considered equilibrated and break-up fragments are built-up from the detected ones attributing ran- 
domly neutrons and light charge particles to heavier fragments {Z > 5) according to pre-defined 
evaporated particle multiplicities and weighted by final sizes of fragments. For the percentages 
of evaporated deuterons, tritons, ^He and ^He we used the following values: pd=0.36, pt=OAQ, 
p4He=0.60 and p6He=0.30 obtained using particle-fragment correlation measurements IJ] and cor- 
recting for pre-equilibrium using the method proposed in Ref. ^. The percentage of evaporated 
protons Pp is set to 0.30 assuming that the contribution of fragments with Z > 35, not consid- 
ered in Ref. [4], modifies exclusively the proton yield. The percentages of evaporated neutrons 
and nuclei with Z = 3, 4 are free parameters which have to be determined simultaneously with 
the average excitation energy of the primary products such as to recover, after sequential particle 
evaporation, the experimental charge distribution. The apart case of ^He whose average kinetic 
energy exceeds by 15 MeV that of tritons and by 13 MeV that of '^He and the urge of not intro- 
ducing extra hypothesis in addition to the break-up scenario, make us tentatively assume that ^He 
originates prevalently from the break-up stage (p3He=0) {z]. Though not relevant for reconstruction 
of break-up fragment partition due to the relatively low multiplicity of ^He {Y{^He)=0.9), if true, 
this assumption may provide an additional test for the accuracy of break-up kinetic properties 
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determination, as we shall see later on. 

We recall that the paradoxical behavior of ^He is already notorious as evidenced in a multitude 



of nuclear reactions covering a broad energy domain {eg. 7.5 GeV/c p+ C, ' Sn 



Ne+U at 



250 and 400 MeV/nucleon [9], 202 MeV/c p+^^C, '^^Ca,^^Cu, ^^'SSMoj^^^U [Iflj], central collisions of 
Au+Au at 100, 150 and 250 MeV/nucleon ll| and Xe+Sn at 50 MeV/nucleon and impossible 
to understand within a thermal scenario (according to which < K >th (^He) ~< K >th (^He)) 
eventually coupled with a collective expansion (providing < K >coii (^He) << K >coii (^He)). The 
high kinetic energy of '^He was so far explained as a signature of its early synthesis from coalescing 
non-equilibrium nucleons |l2] or a time expanding source y^]. As the coalescence as a mechanism 
for fragment formation is expected to work for energies higher than the one which characterizes 
our reaction and the fragment emission from an expanding source [l^ is incompatible with the 
adopted break-up scenario, the assumption of the break-up genesis of ^He, in contradistinction to 
other charged particles and fragments which may originate from both break-up and evaporation 
stages, seems to be the most appropriate. 

So, assuming for the primary fragments a roughly constant excitation energy per mass unit (e) 
y, [l^ , we allow these fragments to de-excite following a standard Weisskopf scheme where realistic 
evaporation Coulomb barriers and complete experimental information on excited states of 

nuclei with ^ < 14 have been implemented. The value of e which reproduces the experimental 
charge multiplicity within 10% is 2.5 MeV/nucleon, in good agreement with the value proposed 
by Hudan et al., 2.2 MeV/nucleon [^]. The "reconstructed" break-up and asymptotic charge 
distributions are presented in Fig. [T]in comparison with the experimental data. For the percentage 
of evaporated neutrons and nuclei with Z = 3,4 we obtained 0.40, 0.14 and 0.10, respectively. 

To complete the static picture of the break-up stage we mention that the average values of total 
multiplicity (Ntot), charged product multiplicity (Nc) and fragment multiplicity (A^z>5) are 30.5 
(asymptotic 53.8), 13.8 (asymptotic 23.2; experimental 20.3) and 4.2 (asymptotic 4.0; experimental 
4.2) respectively. The average total excitation energy {Eex = Yli ^i) 396.9 MeV and Q-value 
equals 126.7 MeV. Both total multiplicity and Q-value are subject to variation if final fragments 
isotopic composition differs from the presently assumed one. 
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FIG. 1: Reconstructed break-up and asymptotic charge distributions corresponding to the (198,83) 
equihbrated source formed in Xe+Sn at 32 MeV/nucleon multifragmentation reaction in compari- 
son with experimental data. 

III. ATTEMPTS TO DETERMINE THE FREEZE-OUT VOLUME AND KI- 
NETIC PROPERTIES OF THE BREAK-UP CONFIGURATION 



The restoration of the fragment configuration event by event suggests the possibihty to deter- 
mine the average values of the other observables which characterize the break-up (the freeze-out 
volume and the thermal and collective energies) by comparing the simulated asymptotic kinetic 
energy distributions < K > vs. Z and dN/dK vs. K with the experimental ones. The idea of 



extracting information on the freeze-out stage and in particular on the freeze-out yo 
energy gain during Coulomb propagation was advanced for the first time in Refs 



181, 



ume using the 
under 



IS 



the simplifying assumptions of low excitation of primary fragments and negligible radial collective 
flow. To reach our goal, we assume full decoupling between fragment formation at freeze-out and 
amounts of energy dissipated in different degrees of freedom and define the freeze-out volume as 
the spherical container in which the primary fragments are localized. Interesting enough, present 
consideration of volume agrees equally with the two extreme scenarios in the literature: volume 
as an external constraint as routinely employed by microcanonical models 2l[ or simply the spa- 



tial extension of the expanding system when fragments cease to interact by nuclear force or final 



multiplicity is attained as in dynamical models 



221. 



For simplicity, fragments are treated as normal nuclear density non-overlapping spheres, their 
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size being thus the experimentally measured one if Z < 6 [24] or calculated according to 
R = 1.2A^/^ for heavier nuclei and the freeze-out volume is assimilated with the smallest sphere 
which includes entirely the fragments. Once localized, fragments share the available thermal ki- 
netic energy, a free parameter in our simulation, according to a Maxwell distribution under total 
momentum and angular momentum conservation laws and then propagate in the free space un- 
der the action of collective motion and Coulomb repulsion. Several values of the freeze-out volume 
have been considered, ranging from 4Vq - the lowest value allowed by the fragment non-overlapping 
condition up to 12Vq - slightly exceeding the maximum evaluations in the literature 



22 



23( 1 ■ where 



Vq stands for the volume at normal nuclear density. Taking into account that recent experimental 
works 



251 ] indicate for the flow profile values larger than 1, we consider this quantity as a free 
parameter. In each case average total thermal energy Kth, magnitude Efiow and profile afiow of 
collective flow have been tuned to provide the best agreement between the simulated asymptotic 
< > vs. Z and dN/dK vs. K distributions and the corresponding experimental data. The 
results indicate that Coulomb and flow energy, both dependent on fragment size and position, may 
compensate each other so that equally good results are produced by any suitably tuned (freeze-out 
volume, collective energy) set. This result is visible in Figs. [2] and [3] where the asymptotic < K > 
vs. Z and dN/dK vs. K distributions of light charged fragments are plotted for three situa- 
tions: (i) V = 4Vb, Efiow=0.68 MeV/nucleon, a//o^„=1.3; (ii) V = GVq, Efiow=0.9 MeV/nucleon, 
a/«ou)=l-2 and (iii) V = 12Vo, -E'/ioii;=l-24 MeV/nucleon, afiow=^-'^- As one may notice from Fig. 
m the agreement between the asymptotic stage simulated < K > vs. Z distribution and the exper- 
imental data is good up to Z=40, while serious underestimation is produced for heavier fragments. 
The reason may be related to the arbitrary way in which break-up fragments larger than Z=35 
have been built-up from the experimental partition in lack of any estimation of the evaporated 
particle multiplicities. In what regards the dN/dK vs. K distributions, in general the simulated 
spectra of light charged particles describe well the experimental ones (see Fig. [3]), while systematic 
underestimation (of about 70 %(for V = 4Vo) - 80 % (for V = 12Vo)) of the observed widths is 
realized for fragments (not shown). 

The only kinetic quantity whose univocal determination was possible following the above pro- 
cedure is the total thermal energy as it influences strongly the slope of the decreasing tail of the 
kinetic energy spectra of light charged particles which is practically insensitive to Coulomb repul- 
sion and collective motion in the considered domains. The obtained value is Xt/i=396 MeV and 
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FIG. 2: Asymptotic simulated average kinetic energy distributions as a function of fragment 
charge corresponding to different values of freeze-out volume and collective energy in comparison 
with experimental data of the Xe+Sn at 32 MeV/nucleon multifragmentation reaction. Complete 
particle evaporation was performed after full propagation, at 500 fm/c. 

corresponds to a kinetic temperature, T = Kth/{3/2Ntot), of 8.65 MeV. We stress here that the 
occurence of nearly identical values for Kth and E^x is fortuitous and should not be assigned a 
physical meaning. In what regards the dN/dK vs. K distributions of fragments, the results of our 
simulation indicate significant and complex dependence on all input quantities (i^th, -E/zoto) ^)) a 
quantification of their effects being impossible. 

It is interesting to notice here that the simultaneous good description of the kinetic energy 
spectra of ^He and of those of p, d, t and ^He supports the hypotheses of (prevalent) break-up 
formation of this nucleus and, at the same time, may be regarded as a validity test of the procedure 
employed to determine Kth- 

A criticism to our procedure may concern the stability of the solution with respect to the 
moment when excited fragments suffer sequential evaporation. To clarify this issue kinetic energy 
spectra have been compared in two extreme cases: (1) full particle evaporation occurs after complete 
propagation, at 500 fm/c and (2) full particle evaporation takes place just after break-up, at fm/c. 
In the second case, the final fragments are propagated in the correspondingly altered Coulomb field 
during 500 fm/c. In particular, the relatively advanced fragmentation of the system at break-up 
and the rather poor particle production during sequential evaporation determine an increase of the 
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FIG. 3: Asymptotic simulated kinetic energy spectra of light charged particles corresponding to 
different values of freeze-out volume and collective energy in comparison with experimental data 
of the Xe+Sn at 32 MeV/nucleon multifragmentation reaction. Complete particle evaporation was 
performed after full propagation, at 500 fm/c. 

total Coulomb energy by only 1 % in the case of early evaporation with respect to the break-up 
stage. The most important effect was evidenced for protons and the dN/dK vs. K distributions 
are plotted in Fig. |4] for the particular case of = 6Vo, -E'coiJ=0.9 MeV/nucleon and acoU=^-'^- 
As one may see in the third panel, the two distributions differ only in the low energy domain, 
the decreasing tail used to determine Kth being practically unmodified. The case corresponding 
to early de-excitation gives a better reproduction of the experimental data and an average kinetic 
energy of Z = 1 particles exceeding by 1.6 MeV the result predicted by the opposite case. The 
influence of the evaporation instant on the kinetic energy spectra may be understood taking into 
account the behavior of the break-up and evaporated proton spectra. Thus, the already discussed 
poor sensitivity of the break-up protons kinetic energy to moderate variations of the Coulomb field 
leads to perfectly superimposable dN/dK vs. K distributions corresponding to the two evaporation 
scenarios. For the sake of completeness, the left panel of Fig. [J] depicts also the dN/dK vs. K 
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FIG. 4: Simulated kinetic energy spectra corresponding to V = GVq, EcoU=0.9 MeV/nucleon and 
acoii = ^-2 for the break-up (left), evaporated (middle) and break-up + evaporated protons (right) 
when sequential evaporation occurs at 500 fm/c (1) and, respectively, fm/c (2). In the third panel 
the experimental distribution is plotted with histogram. In the first panel the final distributions 
corresponding to the two evaporation instances (dotted and dashed lines) nearly coincide pointing 
out the modest increase of the Coulomb energy caused by an early evaporation which practically 
takes place in the freeze-out volume. 

distribution of protons at break-up (t=0 fm/c). Apart the obvious increase of the kinetic energy, 
the propagation in the Coulomb field produces a partial suppression of the low energy range. In 
what regards the evaporated protons, the middle panel of Fig. U] indicates a significant increase of 
the kinetic energy in the case of early evaporation with respect to the (1) scenario. However, the 
widths of these distributions determined mainly by the break-up fragments excitation energy are 
small enough so that only the low energy range of the total dN/dK vs. K distributions proves 
sensitive to the evaporation time (right panel of Fig. 

Once estimated the average energy dissipated in different degrees of freedom, one may evaluate 
the total available energy of the source. Thus, the calorimetric equation, 

Etot = Q + Eex + Kth + y Coulomb + Efiow, (1) 

predicts for all considered situations a total energy of 7.4 MeV/nucleon, the E*={Q + E^x + Kth + 
Vcouiomb) part ranging from 6.15 MeV/nucleon at 12Vo to 6.7 MeV/nucleon at AVq. 
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IV. ADDITIONAL INFORMATION ON THE BREAK-UP STAGE: COMPARI- 
SON WITH A MULTIFRAGMENTATION MODEL 



The limitations of the above procedure to completely characterize the break-up stage of a nu- 
cleus using all available experimental data should not be attributed to the technical flaws of the 
adopted recipe but rather reflect the simplified treatment of the phenomenon and require further 
consideration. Thus, the break-up is defined as the ultimate stage of the cooling-expansion process 
suffered by the hot and initially compressed matter formed in heavy ion collisions when it ceases 
to exist as a mononuclear configuration and crashes into pieces. With this picture in mind, it 
is natural to admit that all break-up observables should depend on the initial energy transferred 
into the system. Recent results obtained within a microscopic approach 20] are interesting in this 



sense as they show in the pre-break-up stage a monotonic increase of both nuclear volume and 
expansion energy with the excitation energy. On the other hand, statistical models which compute 
the weights of all possible fragment configurations C : {{Ai, Zi^ Cj, rj, pj}, i = 1, NF} compatible 
with a certain state of the equilibrated source (eg. for microcanonical models, {A, Z, E, V)) show 
a strong interplay between all energetic degrees of freedom, including Q. In this perspective, one 
possibility to select the unique physical break-up state from the multiple solutions obtained by 
the restoration algorithm is to perform a consistency test with predictions of a multifragmentation 
model, statistical or dynamical. Fragments treatment as hard non-overlapping spheres charac- 
terized by a normal-density zero-temperature binding energy which may be additionally excited 
together with a freeze-out volume configured as a spherical container in which break-up fragments 
are localized recall the break-up stage description of the MMM version 23] of the microcanonical 



multifragmentation models and encourage comparison with the predictions of this model. We re- 
call that, nevertheless, within the microcanonical models, the freeze-out volume has more dramatic 
consequences as it affects not only the kinetic energy spectra, but the fragment partitions them- 
selves via geometric restrictions and, again. Coulomb energy which this time enters the statistical 
weight of a configuration via energy balance |23f]. Finally, as only mean values will be considered 
along this paper, it is reasonable to assume that the freeze-out volume is fixed, without refuting the 
fact that in experiment it may in principle vary from one event to another. A reasonable agreement 
may constitute an additional proof in favor of the break-up scenario. 

Therefore, we start from the same equilibrated source (^0)-^o)=(198,83) and consider for the 
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freeze-out volume and excitation energy the values resulted from the break-up reconstruction al- 
gorithm: {V = Wo, E*=6.7 MeV/nucleon), {V = 6Vb, E*=6.5 MeV/nucleon) and {V = UVq, 
E*=6.15 MeV/nucleon). Because for micro canonical models the break-up stage characteristics 
were proved to depend dramatically on primary fragment definition in terms of excitation energy 
[isl . two situations will be addressed. In both cases fragments are allowed to absorb excitation up 
to the binding level, but in the first case the high energy decrease of the nuclear level density, 



/5(e) = ^^^1/4^5/4 exp(2Vae) exp(-e/r), (2) 
with a = 0.114A + 0.098A^/^ MeV~^ is modulated by a constant limiting temperature factor (r=9 



MeV) 



27[, while a mass dependent one is used in the second case {t{A) = 4.3 + 18.3/\/]4 — 



19.9/^ MeV) 



29(] . Although finite values of r result in different values of internal and kinetic 



temperature, this does not represent an equilibrium violation [30]. The presence of a high energy 
cut-off on the level density is required by the fragments binding criterion and is a common feature 
of models which use clusters and not particles as degrees of freedom. 

Table I summarizes the values of different partial energies and product multiplicities at break- 
up obtained in each considered case, while Fig. [5] depicts the charge multiplicity distributions 
produced by MMM against the reconstructed experimental break-up distribution. 

As one may notice, contrary to what we have seen above, microcanonical models prove extremely 
sensitive to the freeze-out volume - excitation energy (E*) correlated variation. By favoring asym- 
metric configurations with lower fragment multiplicity, a decreasing freeze-out volume leads to a 
broadening of the charge distribution and the drop of the Q-value. Quite interesting is the fact 
that while the evolution of product multiplicities and partial energies is similar for the two consid- 
ered hypotheses of internal excitation, the effect on Y{Z) distributions differs significantly. Thus, 
by increasing the freeze-out volume, a constant limiting temperature produces a reduction of the 
charge distribution range by 20 units (open circles), while for a mass dependent r the modifica- 
tion is ten times smaller (open diamonds) [l^. In what regards the variation of total Coulomb 
repulsion, the values in Table I indicate that the effect of a more uniform population accomplished 
for large volumes due to an increased fragmentation is annihilated by the infiuence of the volume 
size such that an increasing volume produces a monotonic decrease of Vcouiomb- Finally, according 
to MMM more expanded configurations are characterized by relatively colder fragments (e=1.5 - 
2 MeV/nucleon at 12Vo and e=2.2 - 3.2 MeV/nucleon at 4Vo) with correspondingly lower total 
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FIG. 5: Break-up stage fragment charge multiplicity distributions obtained by MMM for different 
states of the source (^0;'^o)=(198,83) as indicated on the figure in comparison with the experi- 
mental reconstructed distribution corresponding to Xe+Sn at 32 MeV/nucleon multifragmentation 
reaction. 

internal excitation. As in microcanonical models thermal kinetic energy is determined event by 
event by the available energy, the variation of Kth results from the interplay between the other 
energetic degrees of freedom and, as indicated in table I, it manifests a monotonic decrease with 
the volume. 

Turning back to the possibility to identify the break-up stage performing a consistency test 
with the predictions of MMM, one may say that even if for V = 4Vo, £'*=6.7 MeV/nucleon 
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TABLE I: Break-up stage values of average total thermal kinetic energy, excitation energy. Coulomb 
energy, Q value and excitation energy per mass unit of primary fragments expressed in McV and 
total multiplicity, charged product multiplicity and fragment multiplicity obtained by MMM for 
the multifragmentation of the (198, 83) source with different freeze-out volumes and excitation 
energies and two definition of primary fragments in terms of excitation energy, as indicated in the 
first column in comparison with the reconstructed experimental break-up stage values. 



(F/Vb, E* (MeV),T (MeV)) 




^^oulomb 


Q 


e/A Ntot N^h Nz>5 


(12, 6.15, 9) 


358.1 269.8 


295.7 


294.0 


2 


37.9 23.4 


4.8 


(12, 6.15, t{A)) 


369.9 196.6 


310.4 


340.7 


1.5 


40.0 25.4 


5.4 


(6, 6.5, 9) 


321.7 425.8 


358.8 


180.6 


2.7 


26.6 17.5 


5.1 


(6, 6.5, t{A)) 


346.5 307.6 


388.0 


244.8 


1.8 


29.2 20.0 


6.2 


(4, 6.7, 9) 


282.4 565.6 


382.2 


96.3 


3.2 


19.5 13.4 


4.8 


(4, 6.7, t{A)) 


316.0 408.3 


435.4 


166.8 


2.2 


22.1 16.0 


6.6 


(4, 6.7, 6.8) 


371.3 401.1 


394.7 


159.5 


2.5 


23.6 16.5 


4.7 


reconstr. exp. break-up 


396.0 396.9 




126.7 


2.5 


30.5 13.8 


4.2 



and r=9 MeV the calculated Y{Z) distribution agrees well with the reconstructed experimental 
one, the result is unacceptable due to the considerably different sharing of total energy: Kfh and 
Q are underestimated by 30% and, respectively, 24 % and E^x is overestimated by 40 %. An 
extra diminishing of the freeze-out volume can not be considered a solution as it will lead to Y{Z) 
distributions broader than the experimental one and Kth, Eg^ and Q will be still worse. The amount 
of energy dissipated in Kth, Eg^ and Q in better agreement with the estimations of the break-up 
reconstruction algorithm obtained for V = 4Vb and E*=Q.7 McV/nucleon when t{A) ranges from 
5.9 (for ^=100) to 8.5 (for ^=5) suggests as a possible solution of our problem the reduction of the 
constant value of the limiting temperature. Thus, diminishing r to the arbitrary value 6.8 MeV, 
we get an Y{Z) distribution which follows closely the "experimental" break-up curve and partial 
energies in acceptable agreement with the "measured" ones. The only quantity for which MMM 
and the break-up reconstruction algorithm give different results is the total multiplicity and the 
explanation relies on the different isotopic composition of the primary fragments in the two cases. 
Indeed, in MMM the break-up products are rather neutron rich and relatively few neutrons are 
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created at break-up. The good reproduction of kinetic distributions < K > vs. Z and dN/dK vs. 
K after sequential evaporation does not require further consideration as it is guaranteed by the 
procedure adopted to build the break-up stage from experimental data. 



V. DISCUSSION AND CONCLUSIONS 

Even if promising as an approximate determination of an average freeze-out volume consistent 
with some relevant experimental information, present results depend upon fragment definition, 
further studies on the shape, density and related internal excitation of nuclei at break-up being 
mandatory for a definite answer. For instance, the ambitious task of checking the consistency 
between statistical and dynamical model predictions is expected to provide information on frag- 
ment characteristics at break-up and answer the fundamental question whether full equilibration 
is attained in heavy ion reactions. A remarkable attempt in this sense, though restricted to frag- 
ment partitions and average kinetic energy as a function of charge and dealing with a slightly 
different source, has been recently published in Ref. [sil reaching the conclusion that a distinct 
equilibrated stage may be identified in the dynamical evolution of the presently considered Xe+Sn 
at 32 MeV/nucleon reaction as described by the Stochastic Mean Field approach [s^]. 

The difficulty to restore the break-up stage using exclusively asymptotic experimental informa- 
tion and the important role played by fragments internal excitation are well illustrated by Refs. 

0,0]. 



For instance, Ref. [15| provides an example on how the agreement between experimental 
average properties of fragments and asymptotic stage predictions of a statistical multifragmenta- 
tion model as criterion for identifying the equilibrated state of a nuclear source may provide for the 
freeze-out volume values differing by a factor of two when various parameterizations of the break- 
up fragments limiting temperature are employed. Even more eloquent is the case of Ref. where 
the assumed relation between the thermal kinetic temperature, break-up fragment internal tem- 
perature and limiting temperature {1/Tfrag = 3/2 < Kjo +1/''") [33] fixes the sharing between 



fragment excitation and thermal energy. This fact, coupled to some differences in the building of 
the freeze-out stage such as the distance between fragments' surfaces, which in the present study 
may approach zero, determines the variance between break-up fragments characteristics obtained 
in Ref. with respect to the present work. 

In conclusion, for Xe+Sn at 32 MeV/nucleon multifragmentation reaction fragment partitions at 
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break-up were built event by event from the experimental fragment partitions using experimentally 
evaluated average multiplicities of evaporated light charged particles which also allow implicit 
estimation of the limited average excitation energy of primary fragments. Kinetic characteristics 
of the break-up stage were inferred from the comparison of simulated asymptotic < K > vs. 
Z and dN/dK vs. K distributions with the experimental ones. Kth was determined uniquely 
from the kinetic energy spectra of light charged particles. Very interestingly, the hypothesis on 
prevalent break-up genesis of ^He adopted to explain the high value of the average kinetic energy 
is supported by the ability to reproduce its broad dN/dK distribution and, at the same time, 
validate the procedure used for the break-up reconstruction. We mention here that the relatively 
low excitation of the source and the numerous proofs in favor of existence of the break-up stage make 
earlier scenarios on ^He formation from coalescing non-equilibrium nucleons or an expanding source 
less plausible. The size of the average freeze-out volume and the corresponding amount of collective 
flow may be selected from the multiple solutions provided by the break-up restoration algorithm 
performing a consistency test with a multifragmentation model with identical fragment definition. 
Acceptable agreement between the MMM results and predictions of the reconstruction algorithm 
was obtained for V = 4Vo assuming for primary fragments a constant limiting temperature t=6.8 
MeV which leads to an average excitation energy per mass unit e/A=2.h MeV, in agreement with 
the experimental evaluations of Ref. 4|. 



[1] F. Gulminelli, Ph. Chomaz, O. Juillet, M. J. Ison and C. O. Dorso, in Proceedings for VI 
Latin American Symposium on Nuclear Physics and Applications, Iguazu, Argentina (2005), 
to be pubhshed in Acta Phys. Hung. A ( nucl-th/0511012[ ). 

[2] J. D. Frankland et al, Nucl. Phys. A689, 905 (2001). 

[3] J. D. Frankland et al, Nucl. Phys. A689, 940 (2001). 

[4] S. Hudan et al, Phys. Rev. C 67, 064613 (2003). 

[5] S. Piantelh et al, Phys. Lett. B627, 18 (2005). 

[6] R. J. Charity, Phys. Rev. C 58, 1073 (1998). 

[7] N. Le Neindre, These de doctorat, Universite de Caen (1999), 
|http:/ /tel.ccsd. cnrs.fr/tel-b6003741 \ S. Hudan, These de doctorat, Universite de Caen 
(2001). 



15 



[8] V. I. Bogatin et al, Yad. Fiz. 32, 1363 (1980). 

[9] H. H. Gutbrod et al, Phys. Rev. Lett. 37, 667 (1976). 
[10] W. Markiel et al, Nucl. Phys. A485, 445 (1988). 
[11] G. Poggi et al, Nucl. Phys. A586, 755 (1995). 
[12] W. Neubert and A. S. Botvina, Eur. Phys. J. A 7, 101 (2000). 

[13] R. Bougault et al, in Proceedings of the XXVII International Workshop on Gross Properties 
of Nuclei and Nuclear Excitations, Hirschegg, 1999, edited by H. Feldmeier, J. Knoh, W. 
Noerenberg and J. Wambach (GSI Darmstadt, 1999), p. 24. 

[14] W. A. Friedman, Phys. Rev. C 42, 667 (1990). 

[15] Ad. R. Raduta, E. Bonnet, B. Borderie, N. Le Neindre and M. F. Rivet, Phys. Rev. C 

72,057603 (2005). 
[16] M. F. Rivet et al, Phys. Rev. C 25, 2430 (1982). 
[17] L. Vaz et al, Z. Phys. A 311, 89 (1983). 

[18] Bao-An Li, D. H. E. Gross, V. Lips and H. Oeschler, Phys. Lett. B335, 1 (1994). 

[19] T. C. Sangster, M. Begemann-Blaich, Th. Blaich, H. C. Britt, L. F. Hansen, M. N. Namboodiri 

and G. Peilert, Phys. Rev. C 51 (1995). 
[20] M. Begemann-Blaich et al, Phys. Rev. C 58, 1639 (1998). 

[21] D. H. E. Gross, Rep. Progr. Phys. 53, 605 (1990); J. P. Bondorf, A. S. Botvina, A. S. lljinov, 

I. N. Mishustin and K. Sneppen, Phys. Rep. 257, 133 (1995); S. E. Koonin and J. Randrup, 

Nucl. Phys. A 474, 173 (1987). 
[22] M. Parlog, G. Tabacaru, J.P. Wieleczko, J.D. Frankland, B. Borderie, A. Chbihi, M. Colonna, 

M.F. Rivet, Eur. Phys. J. A 25, 223 (2005) and references therein. 
[23] Al. H. Raduta and Ad. R. Raduta, Phys. Rev. C 55, 1344 (1997); ibid., Phys. Rev. C 65, 

054610 (2002). 

[24] Landolt-Bornstein, A. M. Hellwege , K. H. Hellwege (Eds.), Numerical data and functional 

relationships in science and technology, vol. I, Springer- Verlag, Berlin, 1961. 
[25] A. Le Fevre et al, Nucl. Phys. A735, 219 (2004). 

[26] J. N. De, S. K. Samaddar, X. Vinas and M. Centelles, Phys. Lett. B638, 160 (2006). 
[27] A. S. lljinov, et al, Nucl. Phys. A543, 517 (1992). 

[28] J. B. Natowitz, R. Wada, K. Hagel, T. Keutgen, M. Murray, A. Makeev, L. Qin, P. Smith, 
and C. Hamilton, Phys. Rev. C 65, 034618 (2002) and references therein. 

16 



[29] M. Baldo, L. S. Ferreira and O. E. Nicotra, Phys. Rev. C 69, 034321 (2004). 

[30] Al. H. Raduta and Ad. R. Raduta, Nucl. Phys. A703, 876 (2002). 

[31] A. H. Raduta, M. Colonna, V. Baran, and M. Di Toro, Phys. Rev. C 74, 034604 (2006). 

[32] M. Colonna et al, Nucl. Phys. A642, 449 (1998). 

[33] S. E. Koonin and J. Randrup, Nucl. Phys. A474, 173 (1987). 



17 



